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1 nt redact i on 


This is the finul report for NASA contract No. NASS- J 1456. The 
contract period began on 15 July 1973 and ended IS January 1974. A 
three-month, no-cost extension of this contract ended on 15 April 1974. 

During the contracted time we developed the theory and computer 
programs, based on rigorous electromagnetic theory, required for the an- 
alysis and design of echelle gratings. These gratings are designed for 
instruments that operate in the ultraviolet portion of the spectrum. 

The theory has been developed so that the resulting computer programs 
will be able to analyze deep (up to 30 wavelengths) gratings by includ- 
ing as many as 100 real or homogeneous diffraction orders. The program 
calculates the complex amplitude coefficient for each of the diffracted 
orders. Also, a check on the numerical method used to solve the in- 
tegral equations is provided by a conservation of energy calculation. 

Theory 

For our work we adopted the rigorous vector diffraction theory 
first described by Zaki^ and Kalhor.^ The essentials of 
this theory are contained in our second-year final report dated July 
1973.^ Below we describe the equations resulting from this theory and 
how they are used to obtain the desired amplitude coefficients. Some of 
the equations are reproduced here for ease of reference. 

For the E-parallel state of polarization the incident electric 
field vector is parallel to the grooves of the reflection grating. Ap- 
plication of the boundary conditions for this polarization state yields 
the following equations: 
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where the integrals are to be evaluated on the grating surface so that 
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y = f(x) and y' = f(x')i where f(x) is the function that describes the 
grating profile. Also, 

Ei is the incident field 

Ei +Ejj is the total electric field (incident + diffracted) 

3E/3n is the derivative of the total electric field with 
respect to the grating surface normal 

Am Is the desired amplitude coefficient for the diffrac- 
tion order 

G' (x,y;x',y is the Green’s function of the problem and in 
the plane wave representation assumes the form 


G(x,y;x',y’) 


-ikx sin6i “ -i2irm(x-x')/a +ik cos9m|y-y' 

^ k^os0m 

C3) 


The various diffraction angles, for a given set of operating conditions, 
are given by the well known grating equation 

sin0^ + 

where a is the grating period and A is the wavelength. 

In Eq. (1) the quantity 9H/3n describes the induced current dis- 
tribution on the surface of the grating and is the only unknown. There- 
fore Eq. (1) can be solved for 9E/3n and the answer can be inserted into 
Eq. (2) so that the unknown amplitude coefficients A^ can be determined. 
Below we briefly describe the method that we have chosen to perform these 
calculations. 

We first note that in order to evaluate the continuous functions 
of Eqs. (1) and (2) they must be cast into a discrete form suitable for 
digital computation. This can be done by first noticing that the inte- 
grals in Eqs. (1) and (2) are over a finite interval (0,a) with respect 
to the variable of integration x. A discrete representation of the quan- 
tities appearing in Eqs. (1) and (2) can be obtained by expanding each 
function in terms of its complex Fourier series expansion with a period 
of a. In this manner we may define the following: 
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j=-« 


( 4 ) 
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G-(x,£(x);x-,f(x')) ^ I I Kj,„ e-^2.nx/a 

j = -oo n=-a» 


{1 +[£'(x)]2}*i 

3n 


I Sn 


+i2Trnx/a 
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^-iBinx ^+iYmf(x) 
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( 7 ) 


n=-«> 


and we have used the definitions 


Sm = 

C2n/A) 

sine„ 

Ym = 

(21T/X) 

cos6jn. 


( 8 ) 


The complex-valued Fourier expansion coefficients hj , Sn> ®m n 

are calculated by making use of the orthogonality of the complex expo- 
nentials appearing on the right-hand sides of Eqs . (4) through (8). 

Their detennination in this way assures the faithful representation of 
the continuous functions in the interval (0,a). Substitution of Eqs. (4) 
through (8) into Eqs. (1) and (2) casts the continuous integral equations 
into a set of equivalent matrix equations 

00 

hj = a 'I ,n Sn (9) 

n=-» 

1 ^ 

^ ” 2ym ^ ®m,n Sn* (10) 

n=-~ 

Equations (9) and (10) each represent an infinite set of simultaneous 
equations. The unknown coefficients gj, must be determined from Eq. (9) 
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once the coefficients hj and kj^n been calculated with sufficient 
accuracy, llie values for gn so obtained are then used in Eq. (10) to 
calculate the unknown amplitude coefficients by simple matrix 
multiplication. 

Computer Program 

We now consider the discrete equations (9) and (10). Although 
the sums involved are over all possible positive and negative integers, 
the actual computer calculation must be limited to a finite number of 
integers. The set of integers chosen depends on the operating conditions 
and the particular grating that is being considered. The set of integers 
must include all of the real orders that are present and a sufficient 
number of evanescent orders that reliable results are obtained. 

Despite the mathematical complexity of the analytical calcula- 
tion, Eqs . (9) and (10) for the digital computation are remarkably sim- 
ple. It is this simplicity that lends this program to be more versatile 
than the Berkeley programs. 

We now discuss some of the features that must be considered when 
calculating the required matrix elements. We begin with the matrix K 
composed of elements Each element must be calculated accurately 

and the dimensionality of K must be properly chosen if the corresponding 
values for gn are to be accurate. In order to calculate the elements 
^j,n grating profile f(x) must be specified. A simple profile is 
shown in Fig. 1. 

y 



Fig. 1. Grating profile f(x). 



5 


For this profile the corresponding expression for the matrix elements 
Kj is given by Eq. (11) (next page). The parameters appearing in Eq. 
(11) are defined as follows: 

A = tana 
B = tanB 

qmn^^^ = - AaYm " ^irn 

+ AaYm ' ^irn 

The expressions for same as for qnm^^^ 

qmn^^^, respectively, except that the slope A is replaced by the slope B 
The parameter xi is the normalized apex location of the profile 

xi = xi/a. 


Note that the expression for ii^volves a sum over m. In order to 

obtain an accurate value of Kj it will be necessary to perform this 

infinite sum unless the terms Kjj^(m) converge rapidly. Let M be the max 
imum value of the integer m and define the corresponding value for Kj 
obtained by using this M as 



M • 

= I Kj (m), 

m=-M 


where K- ^.(m) is the mi^ term in the sum. Let c be a small number. It 

j , II 

is possible to find a value for M such that 


■'j.n 


M 


- 


M-1 


(12) 


When this condition is satisfied, the sum over m is truncated, and it is 
assumed that the resulting value for is the acceptable true value. 

In the computer program we have written, we set M = SO and e = 
10“^. This means that m can assume 101 values (including zero). We 
have found that under these conditions a stable solution for K- „ that 
satisfies Eq. (12) is obtained before m assumes the maximum number of 
allowed values. 
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Despite the length of the expression for Kj (Eq. 11) and the 
realization that the calculation must be performed many times for each 
value of j and n, we have found the program to be very efficient. For 
example, the program for the E-parallel polarization requires less than 
3 sec execution time for 10 diffraction orders in most cases. The exe- 
cution times required to analyze coarse echelle gratings will, of 
course, be considerably longer. However, because the program is effi- 
cient, the required execution times will be reduced considerably. 

Currently we are in the final stages of debugging the computer 
program for the E-parallel polarization. The complexity of the program 
has hampered this effort to some extent, but we have demonstrated that 
our method is indeed workable and that the program is much more effi- 
cient, simpler to run, and capable of analyzing four times as many dif- 
fraction orders as any known existing programs. 

Outlook 

The program for the H-parallel state of polarization has not 
been written even though the theory has been derived. We are awaiting 
the opportunity to completely test the E-parallel polarization before 
proceeding to the H-parallel case. The final combination of the two 
programs will allow complete efficiency evaluations as well as polariza- 
tion studies to be continued. 

The special profile shown below in Fig. 2 has unique properties 
when used in autocollimation. 



Fig. 2. Special grating profile. 
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This profile should be evaluated to verify our previous findings, which 
were based on an approximate theory. We were not able to verify these 
findings by using the Berkeley programs. The following results should 
be established: (a) whether the amplitude coefficients vary as a function 
of groove depth, (b) whether all the energy is found in one order and 
zero or negligible energy is present in all the other orders, and Cc) 
the wavelength interval for which the grating displays a high blaze. 
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